Observation of Fragile-to-Strong Dynamic Crossover in Protein Hydration Water 



o 
o 

(N 



Sow-Hsin Chen, 1 * Li Liu, 1 Emiliano Fratini, 2 Piero Baglioni, 2 Antonio Faraone, 3 and Eugene Mamontov 3 
1 Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, 
Cambridge MA 02139 USA 2 Department of Chemistry and CSGI, University of Florence, 
via della Lastruccia 3, 50019 Florence, Italy 3 Department of Material Science and Engineering, University of Maryland, 
College Park, MD 20742 USA and NIST Center for Neutron Research, Gaithersburg, MD 20899-8562 USA 

At low temperatures proteins exist in a glassy state, a state which has no conformational flexibility 
and shows no biological functions. In a hydrated protein, at and above 220 K, this flexibility is 
restored and the protein is able to sample more conformational sub-states, thus becomes biologically 
functional. This 'dynamical' transition of protein is believed to be triggered by its strong coupling 
with the hydration water, which also shows a similar dynamic transition. Here we demonstrate 
experimentally that this sudden switch in dynamic behavior of the hydration water on lysozyme 
occurs precisely at 220 K and can be described as a Fragile-to-Strong dynamic crossover (FSC). At 
FSC, the structure of hydration water makes a transition from predominantly high-density (more 
fluid state) to low-density (less fluid state) forms derived from existence of the second critical point 
at an elevated pressure. 
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Without water, a biological system would not func- 
tion. Dehydrated enzymes are not active, but a single 
layer of water surrounding them restores their activity. 
It has been shown that the enzymatic activity of pro- 
teins depends crucially on the presence of at least a min- 
imum amount of solvent water 0, 0] . It is believed that 
about 0.3 g of water per g of protein is sufficient to 
cover most of the protein surface with one single layer 
of water molecules and to fully activate the protein func- 
tionality. Thus, biological functions [3j, such as enzyme 
catalysis, can only be understood with a precise knowl- 
edge of the behavior of this single layer of water and 
how that water affects conformation and dynamics of 
the protein. The knowledge of the structure and dynam- 
ics of water molecules in the so-called hydration layer 
surrounding proteins is, therefore, of utmost relevance 
to the understanding of the protein functionality. It is 
well documented that at low temperatures proteins ex- 
ist in a glassy state 0, Q , which is a solid- like structure 
without conformational flexibility. As the temperature 
is increased, the atomic motional amplitude increases 
linearly initially, as in a harmonic solid. In hydrated 
proteins, at approximately 220 K, the rate of the am- 
plitude increase suddenly becomes enhanced, signalling 
the onset of additional anharmonic and liquid-like mo- 
tion [ESSES- This 'dynamical' transition of proteins 
is believed to be triggered by their strong coupling with 
the hydration water through the hydrogen bonding. The 
reasoning is derived from the finding that the protein hy- 
dration water shows some kind of dynamic transition at 
the similar temperature |10L lllj| . Here we demonstrate, 
through a high-resolution quasi-elastic neutron scatter- 
ing (QENS) experiment, that this dynamic transition of 
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hydration water on lysozyme protein is in fact the Fragile- 
to-Strong dynamic crossover (FSC) at 220 K, similar to 
that recently observed in confined water in cylindrical 
nanopores of silica materials [l2l Il3j | . Computer simula- 
tions on both bulk water |14| and protein hydration water 
around lysozyme 15] have led to the interpretation of the 
FSC as arising from crossing the locus of maximum in the 
correlation length ( "Widom line" ) which emanates from 
a critical point into the one-phase region; if this interpre- 
tation is correct, then our experiments provide evidence 
supporting the existence of a liquid-liquid critical point 
in protein hydration water, which previously has been 
proposed only for bulk water |16| . 

Water molecules in a protein solution may be classified 
into three categories: (i) the bound internal water, (ii) 
the surface water, i.e. the water molecules that interact 
with the protein surface strongly, and (iii) the bulk water. 
The bound internal water molecules, which occupy inter- 
nal cavities and deep clefts, are extensively involved in 
the protein-solvent H-bonding, and play a structural role 
in the folded protein itself. The surface water, or usually 
called the hydration water, is the first layer of water that 
interacts with the solvent-exposed protein atoms of differ- 
ent chemical character, feels the topology and roughness 
of the protein surface, and exhibits the slow dynamics. 
Finally, water, which is not in direct contact with the pro- 
tein surface but continuously exchanges with the surface 
water, has properties approaching that of bulk water. In 
this letter, we deal with dynamics of the hydration wa- 
ter in a powder of the globular protein lysozyme. This 
hydration water is believed to have an important role in 
controlling the bio-functionality of the protein. 

The biochemical activity of proteins also depends on 
the level of hydration. In lysozyme, enzymatic activity 
remains very low up to a hydration level (h) of 0.2 (h is 
measured in g of water per g of dry protein) and then in- 
creases sharply with an increase in h from 0.2 to 0.5 1171 . 
Various experiments [lsf and computer simulations |19| 



2 



have demonstrated the strong influence of the hydration 
level on protein dynamics. 

It has been found that many proteins exhibit a sharp 
slowing down of their functions (kinetic of biochemical 
reactions) at a temperature somewhere within the in- 
terval of T between 200 K and 250 K. An analysis of 
the mean-squared atomic displacement, (x 2 ), by u sing 
Mssbauer [20|. X-ray 0, and neutron scattering [2jJ 
spectroscopy, in hydrated proteins shows sharp changes 
around a certain sample temperature range: (x 2 ) varies 
approximately linearly as a function of T at low T 
and then increases sharply above T between 200 K and 
250 K. The sharp rise in (x 2 ) was attributed to a certain 
dynamic transition in biopolymers at this temperature 
range. The coincidence of the characteristic tempera- 
tures, below which the biochemical activities slow down, 
and the on-set of the dynamic transition, suggests a di- 
rect relation between these two phenomena. It has also 
been demonstrated that the dynamic transition can be 
suppressed in dry biop olymers, or in biopolymers dis- 
solved in trehalose [231 - It can also be shifted to higher 
temperatures, e.g. between 270 K and 280 K, for proteins 
dissolved in glycerol [2]|. Thus the solvent plays a cru- 
cial role in the dynamic transition in biopolymers. This 
observation led to a suggestion by many authors that 
proteins are 'slaves' to the solvent |24|. Despite many 
experimental studies, the nature of the dynamic transi- 
tion in proteins remains unclear. Many authors interpret 
the dynamic transition as a kind of glass transition in a 
biopolymer 0, S ■ O ur experiments described below 
demonstrate that the origin of the characteristic temper- 
ature controlling both the activity of the protein and the 
transition in the behavior of (x 2 ) is the FSC phenomenon 
in the hydration water, which shares the same crossover 
temperature with the protein. 

Using high-resolution QENS method and the Relaxing- 
Cage Model (RCM, described in Methods section) [2f| 
for the analysis, we determine the temperature depen- 
dence of the average translational relaxation time, (tt), 
for the hydration water. The dynamic crossover tem- 
perature of hydration water is defined as follows. At 
high temperatures, (tt) follows a super- Arrhenius be- 
havior (called a 'fragile' behavior j27| ) describable ap- 
proximately by a Vogel-Fulchcr-Tammann (VFT) law: 
{vr) = Ticxp [DT n /(T — T )], where D is a constant pro- 
viding the measure of fragility and To, the ideal glass 
transition temperature at which the relaxation time ap- 
pears to diverge. In reality, however, this divergence is 
avoided by the system. Instead, an Arrhenius behav- 
ior (called a 'strong' behavior 27]) sets in below the 
crossover temperature Tl where the functional depen- 
dence of the relaxation time switches to a law: (tt) — 
T\exp[EA/ RT]. In this equation, Ea is the activation 
energy for the relaxation process and ks, the Boltz- 
mann constant. The crossover temperature Tl is de- 
fined by the intersection of these two laws, which gives 
1/Tl = X/Tq— (DkB)/EA- In the case of hydration water 
in lysozyme, we found Tl = 220 K, which agrees well with 
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FIG. 1: Measured neutron spectra. Panels A and B show 
normalized QENS spectra at Q = 0.87 A -1 , displaying the 
height of the peak (A) and the wing of the peak (B), respec- 
tively, at a series of temperatures. The inset C shows a plot 
of the peak heights versus temperature. Arrow signs in B are 
intended to highlight the wing of the spectra at the crossover 
temperature. 



the characteristic transition temperature in protein ob- 
served before Since the average relaxation time (ry) 
is a measure of the mobility of a typical hydration water 
molecule, this result implies that the sudden change in 
the trend of mobility of water molecules at the crossover 
temperature triggers the so-called glass transition of pro- 
tein molecules HUEHH. 



Results 

We show in Fig. 2] as an example, a complete set (tem- 
perature series) of QENS area-normalized spectra taken 
at Q — 0.87 A" 1 at ambient pressure. The broaden- 
ing of the quasi-elastic peaks at the wing becomes more 
and more noticeable as temperature increases. And at 
the same time, the peak height decreases accordingly be- 
cause the area is normalized to unity. In the inset C, 
we plot the peak height as a function of temperature. It 
is noticeable that the rate of increase as a function of 
temperature is different across the temperature 220 K. 
From Fig. ^3, we may notice, from wings of these spec- 
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FIG. 2: Neutron spectra and their RCM analyses. Panel A 
displays measured QENS spectra (solid symbols) and their 
RCM analysis results (solid lines) at Q = 0.87 A^ 1 and at 
a series of temperatures. Panel B singles out one particular 
spectrum at T = 230 K and contrasts it with the resolution 
function of the instrument for this Q value (dash line). 



tral lines, that two groups of curves, 270 K > T > 230 K 
and 210 K > T > 160 K, are separated by the curve at 
a temperature 220 K. This visual information obtained 
from the spectra before data analysis reinforces the re- 
sult of the detail line shape analysis to be shown later in 
Fig. |3 that there is an abrupt dynamical transition at 
T L = 220 K. 

Fig-Glshows the result of RCM analyses of the spectra 
taken at Q — 0.87 A -1 for a series of temperatures rang- 
ing from 270 K to 160 K (panel A), and in particular at 
T = 230 K (panel B). In this figure, we display the instru- 
ment resolution function purposely for comparison with 
the measured spectrum. RCM, as one can see, repro- 
duces the experimental spectral line shapes of hydration 
water quite well. The broadening of the experimental 
data over the resolution function leaves enough dynamic 
information to be extracted by RCM. This means that it 
requires a high-resolution backscattering instrument with 
an energy resolution of 0.8 /zeV to adequately study the 
FSC phenomena in hydration water. 

In Fig. |3 we first present (in panel A) the mean- 
squared atomic displacement (a; 2 ) of the hydrogen atoms 
(calculated from the translational Debye- Waller factor, 
Sh(Q,u — 0) = exp[— Q 2 (x 2 )]) versus T to indicate that 
there is a hint of a dynamic transition at a temperature 
between 200 K and 220 K. In panel B, we then present 
in an Arrhenius plot the temperature dependence of the 
average translational relaxation time, (tt), for the hydro- 
gen atom in a water molecule calculated by Equs. (1-2) 
in the Methods section. The contribution from hydro- 
gen atoms in the protein has been subtracted out during 
the signal processing. It is seen that, in the temper- 
ature range 270-230 K, (tt) obeys VFT law, a signa- 
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FIG. 3: Evidence for the dynamic transition. Panel A shows 
the temperature dependence of the mean-squared atomic dis- 
placement of the hydrogen atom at 2 nanosecond time scale 
measured by an elastic scan with resolution of 0.8 fieV. Panel 
B shows temperature dependence of the average translational 
relaxation times plotted in log((Tr)) versus Ta/T, where To 
is the ideal glass transition temperature. In this hgure, there 
is a clear and abrupt transition from a VFT law at high tem- 
peratures to an Arrhenius law at low temperatures, with the 
fitted crossover temperature Tl = 220 K and the activation 
energy Ea = 3.13 kcal/mol extracted from the Arrehenius 
part indicated in the figure. 

ture of fragile liquid, quite closely. But at T = 220 K it 
suddenly switches to an Arrhenius law, a signature of a 
strong liquid. So we have a clear evidence of FSC in a 
cusp form. The T for the fragile liquid turns out to be 
176 K, and the activation energy for the strong liquid, 
E A = 3.13 kcal/mol. 



Discussion and Summary 

Recently, E. Mamontov observed a similar dynamic 
crossover in the surface water on cerium oxide powder 
sample |28j . The surface of cerium oxide is hydroxylated. 
The coverage of water is about 2 layers and the crossover 
temperature is said to be at 215 K. The observed slow 
dynamics is attributed to the effect on the translational 
mobility of the water molecules in contact with the sur- 
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face hydroxyl groups. Thus, our observation of the FSC 
in hydration water of protein may be a universal phe- 
nomenon for surface water. 

It should be noted that the FSC in confined super- 
cooled water is attributed to the crossing of the so-called 
Widom line in the Pressure- Temperature (phase) plane 
in a recent MD simulation work on bulk wate r I14| and 
protein hydration water at ambient pressure 15]. The 
Widom line is originated from the existence of the sec- 
ond critical point of water and is the extension of the 
liquid-liquid coexistence line into the one phase region. 
Therefore, our observation of the FSC at ambient pres- 
sure implies that there may be a liquid-liquid phase tran- 
sition line in the protein hydration water at elevated pres- 
sures. This dynamic crossover, when crossing the Widom 
line, causes the layer of the water surrounding a protein 
to change from the 'more fluid' high-density liquid form 
(which induces the protein to adopt more flexible confor- 
mational sub-states) to the 'less fluid' low-density liquid 
structure (which induces the protein to adopt more rigid 
conformational sub-states). 

In summary, an investigation of the average transla- 
tional relaxation time, or the alpha-relaxation time, of 
protein hydration water as a function of temperature 
reveals a hitherto un-noticed Fragile-to-Strong dynamic 
crossover at 220 K, close to the universal dynamic tran- 
sition temperature documented for proteins in literature. 
This fact implies that the sudden transition of the water 
mobility on the surface of a protein at the FSC triggers 
the so-called glass transition, which is known to have a 
profound consequence on biological function of the pro- 
tein itself. 



Methods 

Sample preparation 

Hen egg white lysozyme used in this experiment was 
obtained from Fluka (L7651, three times crystallized, 
dialysed and lyophilized) and used without further purifi- 
cation. The sample was dried under vacuum in the pres- 
ence of P2O5 to remove any water left. The dried protein 
powder was then hydrated isopiestically at 5 °C by expos- 
ing it to water vapor in a closed chamber until h = 0.3 is 
reached (i.e. 0.3 g H2O per g dry lysozyme). The hydra- 
tion level was determined by thermo-gravimetric analysis 
and also confirmed by directly measuring the weight of 
absorbed water. This hydration level was chosen to have 
almost a monolayer of water covering the protein sur- 
face [2|j. A second sample was then prepared using D2O 
in order to subtract out the incoherent signal from the 
protein hydrogen atoms. Both hydrated samples had the 
same water or heavy water/dry protein molar ratio. Dif- 
ferential scanning calorimetry (DSC) analysis was per- 
formed in order to detect the absence of any feature that 
could be associated with the presence of bulk-like water. 
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FIG. 4: Temperature dependence of the exponents /3 and fly 
giving respectively the power laws oft and Q-dependence of 
ISF. 



Neutron experiments 

High-resolution incoherent QENS spectroscopy 
method is used to determine the temperature depen- 
dence of the average translational relaxation time for 
the hydration water. In this experiment, we measured 
both an H2O hydrated sample and a D2O hydrated 
sample, and take the difference to obtain the signal 
contributed solely from hydration water. Because 
neutrons are predominantly scattered by an incoherent 
process from the hydrogen atoms in water (rather 
than by the coherent scattering process from the 
oxygen atoms), high-resolution QENS technique is an 
appropriate tool for the study of diffusional process of 
water molecules. Using the High-Flux Backscattering 
Spectrometer (HFBS) in NIST Center for Neutron 
Research (NCNR), we are able to measure the average 
translational relaxation time (a-relaxation time) from 
60 ps to 20 ns over the temperature range of 270 K to 
180 K, spanning both below and above the FS crossover 
temperature. For the chosen experimental setup, the 
spectrometer has an energy resolution of ± 0.8 /ieV and 
a dynamic range of ± 11 /ieV [30|. in order to be able 
to extract the broad range of relaxation times covering 
both the fragile and the strong regimes of the relaxation 
times from measured spectra. 

Data Analysis 

QENS experiments essentially provide us with the 
Fourier transform of the Intermediate Scattering Func- 
tion (ISF) of the hydrogen atoms, FniQjt), of water 
molecules in the hydration layer. MD simulations have 
shown that the ISF of both bulk (3l| and confined [3^] 
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supercooled water can be accurately described as a two- 
step relaxation: a short-time Gaussian-like (in-cage vi- 
brational) relaxation followed by a plateau and then a 
long-time (time longer than 1.0 ps) stretched exponen- 
tial relaxation of the cage. The RCM 26], which we 
use for data analysis, models closely this two-step relax- 
ation and has been tested extensively against bulk and 
confined supercooled water through MD and experimen- 
tal data 0, 1^, |32. By considering only the spectra 
with wave vector transfer Q < 1.1 A -1 , we can safely ne- 
glect the contribution from the rotational motion of water 
molecule in ISF j^. The RCM describes the transla- 
tional dynamics of water at supercooled temperature in 
terms of the product of two functions: 



F H (Q,t) 
tt(Q) 



F T (Q,t) = F s {Q,t) exp 
to (0.5QH 



-WMQ)Y 



(1) 



where the first factor, F s (Q,t), represents the short- 
time vibrational dynamics of the water molecule in the 
cage. This function is fairly insensitive to temperature 
variation, and thus can be calculated from MD simula- 
tion. The second factor, the a-relaxation term, contains 
the stretch exponent /3, and the Q-dependent transla- 
tional relaxation time tt(Q), which is a strong function 
of temperature. The latter quantity is further specified 
by two phcnomcnological parameters To and 7, the expo- 
nent controlling the power-law Q-dependence of tt (Q). 
The average translational relaxation time, which is a Q- 
independent quantity, we use in this paper is defined as: 



(r T ) = T r {1/13) //?, 



(2) 



where T is the gamma function. The temperature depen- 
dence of the three phenomcnological parameters, To, /3, 
and 7, are obtained by analyzing simultaneously a group 
of nine quasi-elastic peaks at different Q values. Then 
the average translational relaxation time, (tt), is calcu- 
lated according to Equ. (2) using To and (3. As can be 
seen from Equs. (1) and (2), the product [3j is an expo- 
nent expressing the Q-dependence of the ISF. Fig. Ogives 
the temperature dependence of the exponent (3j, which 
indicates a precipitous drop at the crossover temperature 
220 K from a high temperature value of about 1 to a low 
temperature value of about 0.2, while (3 maintains a value 
of 0.5 all the way through [12|. Note that the values of 
exponents f3 and fij for free diffusion are 1 and 2, respec- 
tively, and the lower values of these two exponents sig- 
nify an anomalous and more restrictive mobility of water 
compared to that of the free diffusion. The Fig. 0] shows 
that while the mobility of hydration water molecule at ps 
to ns timescale deviates significantly from free diffusion 
at all temperatures measured, as temperature goes be- 
low the crossover temperature Tl, the mobility becomes 
drastically reduced. 
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